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Latest Results from the SGO Collaboration 

Presented by S. Collins a * 

a Dept. of Physics and Astronomy, Glasgow University, Glasgow, G12 8QQ, Scotland 

• We present results for the spectrum and decay constants of B mesons from NRQCD using dynamical configu- 

rations at f3 = 5.6 with two flavours of staggered fermions. The light quarks are generated using the Clover action 
with tadpole improvement. 
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1. INTRODUCTION 

We present a study of the heavy-light meson 
spectrum and decay constants, with a focus on 
investigating heavy quark symmetry as well as 
providing quantitative predictions. Our approach 
is to use NRQCD for the heavy quark and in this 
analysis the action includes terms to 0(1/ Mq): 
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Matrix elements must be calculated consistently 
to this order. Here, we restrict the analysis to 
the tree-level operators for the axial and vector 
current to O(l/M ) fl, 
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where T = 7570 and 7, respectively. Tadpole 
improvement is used throughout and then cb is 
given the tree-level value of 1.0. 

The heavy quark propagators were generated 
using the evolution equation 
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for all timeslices. Equation (g) is fully consistent 
to 0(1 /Mo). This differs from the evolution equa- 
tion previously used in B by a factor of 1- 
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the source timeslice, and our results had a system- 
atic error of 0(o,Aqcd/M). However, we found 
this error to be small Q and our previous findings 
are not affected. 

The simulation was performed on 100 16 3 x 32 
lattices at (3 = 5.6 with two flavours of dynam- 
ical staggered fermions, generated by the HEM- 
CGC collaboration. The light quark propagators 
were computed using the Clover action with tad- 
pole improvement, at three masses of light quark, 
hi = 0.1385, 0.1393 and 0.1401. The heavy quark 
propagators were computed over a wide range of 
bare quark mass, 0M0 = 0.8—10.0, in order to 
study heavy quark symmetry. Full details of our 
analysis can be found in 

2. RESULTS 

A vital part of a lattice calculation of B physics 
is a comprehensive calculation of the spectrum. 
This provides a firm foundation for the simu- 
lation of more complicated processes. Figure |l| 
presents the results for the low lying B spectrum; 
a -1 = 2.0 GeV from light spectroscopy has been 
used to convert to physical units and the B me- 
son mass to fix oMq. For the mass splittings 
B s —B and B**—B, the dominant systematic er- 
rors are estimated to be comparable to or less 
than the statistical errors. With the exception of 
B*—B we find good agreement with experiment, 
although, the experimental results for the P—S 
splitting are preliminary. The discrepancy with 
the real world for the hyperfine splitting may be a 
residual quenching effect. However, it is also pos- 
sible our value of cb may be too low by 10—20%; 
a nonperturbative determination of cb is needed. 
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Note that the gross shape of the spectrum 
supports the naive picture of a massive heavy 
quark with velocity u~0.1 surrounded by a 
light quark cloud. In this model the hyper- 
fine spitting is expected to be Mv 2 ~50 MeV, 
while the P— S and 2S—1S splittings are due 
to excitations of the light quark and so are of 
O(A Qcr ,)-200-500 MeV. We can investigate 
this further by quantifying the change in these 
splittings with Mq. Table |l| summarises the co- 
efficients extracted from fitting the spin-averaged 
ground 5-state energy (E) and mass splittings to 
the function 

E, AE = Co + C x /M + C 2 /M 2 + ... (7) 

where M is the PS meson mass. 

The results are compared with our expecta- 
tions in Table |l|. With the exception of the 
bare kinetic energy of the heavy quark (extracted 
from E) the overall agreement is good. While 
(lS|£fcj„ \lS) bare is small, it is opposite in sign 
to the physical result, which is positive. This is 
due to our use of tadpole-improved matrix ele- 
ments Q ; using mean field theory the unimproved 
matrix element corresponds to a large positive re- 
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result is in agreement with an explicit calcula- 
tion by Crisafulli et al || , and demonstrates that 
the matrix element is tadpole dominated. In fact, 
the renormalisation required to obtain the (posi- 
tive) physical quantity from the tadpole-improved 
(IS* I Chin \ ^S)bare 1S smau an d can be calculated 
perturbatively. 

The slope of the 2S—1S splitting gives the dif- 
ference between the physical kinetic energy of 
the heavy quark in the excited state and ground 
state mesons; it is a physical quantity since the 
renormalisation shift cancels in the difference. 
Note that the slope is positive, and its magni- 
tude suggests (2S\ Chin 1 25) is large compared to 
(\S\C km \\S). 

Now consider a similar analysis of the vector 
and axial currents. We define the V and PS de- 
cay constants in the same way using 



fVM = 
S(fVM) = 
(fVM) tot = (fVM) + 5(fVM) 
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Figure 1. The errors shown are purely statistical. 



where 5{f\M) denotes the 0(1/M) correction 
to the current. The expansion of j\f~M about the 
static limit is given in HQET to 0(1 /M) [§ by 



(fVM) tot = {fVM) stat (l + c P /M) (11) 

Cp = Gkin + 2dGhyp + dGcor/6 (12) 

where d — 3 and —1 for PS and V mesons respec- 
tively. In equation ([[2]) the slope is decomposed 
into the contributions from the kinetic energy of 
the heavy quark, Gkin, the hyperfme interaction, 
Ghyp, and the correction to the current, G CO r- 
Equation (JlJ) suggests that these contributions 
can be isolated by finding the slope of various 
combinations of the V and PS decay constants. 

Table || presents our results with the cor- 
responding contributions to the slopes, ob- 
tained using a fit function of the form (0). 

The physical quantities (fy/M)^, (/a/M) and 
(fVilfps/UVMYps are naively expected to 
have coefficients of the slope of O(Aqcd)- How- 
ever, we find a much larger slope than this for the 
total decay constant. The slope of the spin aver- 
age of the decay constants, which only depends on 
Gkim shows that it is the contribution from the 
kinetic energy of the heavy quark which is large 
and dominates cp. Comparing with Neubert Q 
we find good agreement for the slopes of the phys- 
ical combinations of decay constants while for 
the unphysical slope of (/V / M)ps/(/v / M)y and 
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Table 1 

The coefficients of various spectral quantities as determined from first order perturbation theory in l/M, 
where (1S\ represents an 5*-state in the infinite mass limit. The results are in physical units (in powers of 
GeV). Note that those obtained for 2S—1S are estimates and not obtained from a fit. is the energy 
of a heavy quark in the static theory, calculated perturbatively. k; = n c . 



Quantity 


Coefficient 


Results 


Expectation 


E 


C = E°° 


1.00(2) 


Static Result 




A = - E£° 


0.3(3) 






C 1 = (lS\C kin \lS) bare 


-0.1(1) 


+ A 2 


E(2S)-E(1S) 


c = e co (2S)~e co {is) 


0.5 


+^QCD 




d = (25| C km \2S) vhys - (1S\ C km \lS) phvs 


0.3 


+vc 


EfSJ - E('S ) 


C o = 


0.000(2) 







Cl = (1S\ C hyp \lS) phvs 


0.20(4) 


+ A 2 



intercept of 2M S(fVM) PS /(fVM)ps, where 
agreement is not expected, the results differ in 
magnitude and for the former also in sign [Q. 

Extracting the value of the PS decay constant 
at Mb we find 

Z^fB = 216 ±(6) ±(10) ±(20) MeV (13) 

The first error is statistical and the second is 
the systematic error arising from the clover light 
fermions estimated to be 0(a 2 ). The third is due 
to the omission of 0(1/ M^) terms in the NRQCD 
action, and is the dominant error in fg. Using 
the static renormalisation factor of Z^ at = 0.8 
we find / B ~180 MeV. 
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Ci/Co 


Results 




Gkin + dGhyp + G cor /2 


-2.0(4) 


ifVM) 


Ckin 


-1.8(2) 


(fVM)ps 

(fVM) v 


&>Ghyp 


0.72(6) 


(/VM)*?* 
tfVM)*?* 


$>Ghyp + 2G cor /3 


-0.24 


Quantity 


C 


Results 


2M„S(fVM)p S 
(fVM)ps 


G cor 


-1.24(2) 



Table 2 

The results are in GeV. Note that renormalisation 
factors are not included, k; = n c . 



3. CONCLUSIONS 

We have performed a comprehensive study of 
the heavy-light meson spectrum and decay con- 
stants. The question that arises from this analy- 
sis is why there is such a large slope for the decay 
constant when the deviations from the static limit 
for spectral quantities are in agreement with naive 
expectation. A potential model approach || sug- 
gests that this is a consistent picture and that a 
large value for cp is coupled to a large value for 
(2S\C km \2S). 
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